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A variety of physical phenomena can lead to viscous effects. Several sources of shear and bulk 
viscosity are reviewed with an emphasis on the bulk viscosity associated with chiral restoration and 
with chemical non-equilibrium. We show that in a mean-field treatment of the limiting case of a 
second order phase transition, the bulk viscosity peaks in a singularity at the critical point. 

\Q ' I. INTRODUCTION AND THEORY 

o" 

Viscosity has attracted remarkable attention in the RHIC (Relativistic Heavy Ion Collider) community during the 
first years of running. In particular, experimental observations of large elliptic flow have pointed to a small shear 
^"> viscosity and inspired the term "perfect liquid" [lj. In this paper, we review the general theoretical definition of 
viscosity, then show how five different physical effects can lead to non-zero viscous coefficients (Sec. 2). We focus on 
bulk viscosity, and show that one can find large, even singular, effects in the neighborhood of T c (Sec. 3). Although 
the present study focuses on understanding the behavior and physical explanation of the coefficients, we speculate on 
^j ' the experimental manifestations large viscosities might bring about. 

In non-viscous hydrodynamics the elements of the stress-energy tensor depend only on the energy density e and 
£N| particle-number densities ft, 

00 I ^f™W)=PM)<^, (1) 

where P is the pressure and e and n are implicitly functions of r and t. The tilde denotes that Ty is evaluated in a 
frame where the collective velocity u(r) = 0. In Navier-Stokes hydrodynamics, viscosity is incorporated by altering 
T so that it includes terms proportional to the velocity gradients dui/dr 

f i f s - ) =P(e,fr)^+rj(e > Tl)^ + ^-|(V-u)rf,, ) - 1M< . ,7)1 V ■ uM„ .2) 






O 



X 



Here, r\ and B are the shear and bulk viscosities. In ideal hydrodyanmics Vu = (de/dt)/(P+e), so the bulk viscosity 
can be interpreted as describing the correction to the pressure linearly proportional to the rate at which the energy 
density changes, whereas the shear viscosity describes the asymmetry of Ty due to an anisotropic expansion. In non- 
viscous hydrodynamics, accelerations are proportional to the gradient of the pressure, while in general, accelerations 
arise from derivatives of the stress energy tensor, 



c?: ^ +r «>-flf = -ft^- ( 3 ) 

Thus, the components of the stress-energy tensor can be considered as representatives of the pressure in a given 
direction, and any reduction/rise of Ty from viscous effects will result in a slowing/acceleration of the expansion in 
that direction. 

Viscous coefficients can be expressed in terms of correlations in the stress-energy tensor through Kubo relations. 
These are derived by considering alterations of Ta due to a perturbation V . In linear response theory, 

5{T ij (r = 0)) = -(i/n) f d±r'{[T ij (r = 0),V(r')]), V(r') = r^d iUj )T 0j (r'), (4) 

Jr' o <0 

where the perturbation represents the change to the Hamiltonian due to boosting according to a linear velocity 
gradient. 

To derive the Kubo relations, one first makes the substitution Ty — ► ATy, where ATy refers to the difference with 
respect to the time- averaged value of Ty. After inserting (d r 'r' ) into the integral one can integrate by parts, and 
because of the substitution with ATy, one can dispose of the contributions at r' Q — —00 to obtain, 

*<T 4i (r = 0)) = -(i/n) f d 4 r'(d r >r' )([AT ij (r = 0),AT 0j (r')])r' i (d i u j ), (5) 
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= {i/%) / d 4 rVJ([AT tJ (r = 0), ^ AT 0i (r')])rJ(a.« i ), 

= (i/h) [ d i r'r' o ([AT ij (r = 0),AT kl (r')})d k u l . (6) 



The last step involved applying the conservation of the stress-energy tensor, dtTio — —djTij. 

For i 7^ j, symmetries constrain A: and / to equal i and j, which allows the extraction of the shear viscosity from 
Eq. ©, 



r, = (i/h) d*r / r , o<[AT«(0),AT ij -(r , )]) I i^j (7) 

= B,2^/ d4rVU;t ' <[A ^ (0) ' A ^' (r ' )]) - 
By considering the case where 9jUj = (1/3) J^ -V • u, one can inspect I# in Eq. JfjJ) to find the bulk viscosity, 

B=(*/3ft)W d 4 r'r' ([AT rl (0)^T n (r')}) (8) 
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= i™ ^ E / dVe^'dAI^O), AT^(r')]). 

The Kubo relations, Eq.s J7J and JSJl, are fully consistent with quantum mechanics. The classical limit can be obtained 
by first noting that (• • ■) refers to a thermal average with temperature T — 1//3, then applying the identity [2j, 

e-^Vit) = e i/3ndt V(t)e- 0H , (9) 

to one of the terms in the commutator in Eq.s J7J or JSJ, then keeping the lowest term in h, 

Tr e-^AT^O), AT kl (r)} « -ih/3Ti d t AT l0 {Q)AT k i(r), (10) 

which after an integration by parts gives the classical limit of the Kubo relations, 

r, « f d 4 r'(AT y (0)AT i ,(r / )), Mi (H) 
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The classical limit has been applied to determine the shear viscosity from simulations of molecular dynamics 

Although the Kubo relations are difficult to interpret physically, they do make it clear that viscosity is related to 

the size and to the damping of fluctuations of the elements Ty . If fluctuations in Ty (at fixed energy) are large, or if 

they are slow to relax, a large viscosity will ensue. 

Finally, since the Kubo relations are based on linear-response theory, i.e., assuming that the perturbation V(r) is 

small, we emphasize that the corrections to the stress-energy tensor represent an expansion in the velocity gradient. 

If the stress-energy tensor is strongly altered, it brings into question the validity of the linear approximation. This 

would be especially true if the viscous coefficients diverge as will be the case for the example discussed in Sec. IIIII 

II. FIVE SOURCES OF VISCOSITY 

Viscous effects arise whenever the elements of the stress-energy tensor, T^, have difficulty maintaining the equilib- 
rium values due to a dynamically changing system, i.e., one with velocity gradients. In this section we briefly review 
five physical sources of viscosity, the first three of which have already been explained in the literature. Although 
viscosity in non-perturbative systems with ambiguous degrees of freedom might defy the simple descriptions of the 
five effects enumerated below, these examples provide physical insight into the richness of possible root causes for 
shear and bulk viscosity. 



1. Viscosity from non-zero mean free paths: This is the most commonly understood source of viscosity. It 
is straight-forward to see how a non-zero collision time leads to an anisotropy for T^ by considering a velocity 
gradient for a Bjorken expansion, u z — z/t, or equivalently, the velocity gradients are d z u z = 1/r, d x u x — d y u y = 
0. We consider a particle whose momentum is p' z (r) when measured in the frame moving with the collective 
velocity corresponding to its position. In the absence of collisions, p' z will fall with r since the particle will 
asymptotically approach a region where its velocity equals the collective velocity, p' z (t + St) = p' z {t)t/(t + St). 
Meanwhile, p' x and p' are frozen. The resulting anisotropy in the stress-energy tensor yields the following 
expression for the shear viscosity pa, |g , 

r, = (4/5)Pr c , (13) 

where Tr is the collision time. The anisotropy increases the transverse pressure, giving radial flow an initial 
boost |2], and decreases the longitudinal pressure, thus reducing the longitudinal work which results in a larger 
transverse energy p] . It is also easy to see how such an expansion does not yield a bulk viscosity for either ultra- 
relativistic or non-relativistic gases. In those cases an isotropic expansion scales all three momenta proportional 
to 1/r which maintains thermal equilibrium, and collisions do not play a role. This is not the case when m ~ T, 
or especially if the gas has a mixture of relativistic and non-relativistic particles. 

2. Viscosity from non-zero interaction range: If the range of interaction between two-particles extends a 
distance R, interactions will share energy between particles from regions with different collective energies. A 
particle at r = 0, where the collective energy is zero, will share energy with particles whose collective energy is 
(l/2)m(Rd r u) 2 . For Boltzmann calculations, the viscosity will be proportional to PR 2 /t c |g, with the constant 
of proportionality depending on the form of scattering kernel. Both bulk and shear terms result from non-zero 
interaction range. In Boltzmann calculations, the range of the interaction can approach zero for fixed scattering 
rates if the over-sampling ratio is allowed to approach infinity. Although this solves causality problems 9], it 
simultaneously eliminates viscous terms arising from finite-range scattering kernels, which might or might not 
be desirable. This has profound effects on calculations of elliptic flow, which can vary by a factor of 2 depending 
on the range of the scattering kernel [§) . 

3. Classical Electric Fields: Color flux tubes form after the exchange of soft gluons between nucleons passing 
at high energy, and might also be formed during rapid hadronization. Additionally, longitudinal color electric 
fields might be created during the pre-thermalized stage of the collision (color-glass condensate). Since these 
fields tend to align with the velocity gradient, they can be a natural source of shear viscosity. In fact, if the fields 
are purely longitudinal, the elements of the stress-energy tensor become T zz = — e, T xx = T yy = e. Thus, the 
transverse pressure becomes three times as stiff as a massless gas, P — e/3, which is usually considered a stiff 
equation of state. The negative longitudinal pressure signifies that the energy within a given unit of rapidity 
is increasing as the work done by the expansion is negative, similar to the stretching of a rubber band. This 
hyper-shear can lead to the development of early collective radial flow. A sophisticated calculation including the 
effects of interactions amongst the fields 10] showed a somewhat dampened anisotropy compared to the simple 
limit discussed here, with T xx = T yy « 0.5e and T zz rs 0. 

4. Non-equilibrium chemistry: Chemical equilibirum can not be maintained unless the rate at which equilib- 
rium abundances change is much smaller than the chemical equilibration rate, l/r c hemj 

dN 

--Tr- = -(l/T chcm )(N ~ N cq ). (14) 

If the equilbrium number is slowly changing, abundances will vary from equilbrium by an amount, 

SN = -r chcm ^. (15) 

at 

To associate this departure from equilibrium as a viscosity, one must consider the corresponding change in the 
pressure, 
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and make a connection between dN oq /dt in Eq. 1151) with V • u, 
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— r- 1 = -s2— -sv • u. (17) 
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Here, fl is the volume and the second relation exploits the fact that entropy is conserved in a slow expansion. 
The bulk viscosity is then found by comparison of the resulting change in pressure with the definition of viscosity 
inEq. 10), 
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The bulk viscosity will be large whenever the equilibrium number is rapidly changing, e.g., the temperatures are 
falling below the masses, or masses are rising due to restoring chiral symmetry. If the hydrodynamic equations 
explicitly treat particle numbers as current obeying chemical evolution rates, chemical non-equilibration would 
not need to be accounted for through viscous terms. 

5. Viscosity from dynamic mean fields: Bosonic mean fields, such as the a field, obey the Klein-Gordon 
equation. For fluctuations of wave number k — > 0, 

^Aa(t) = -m ff (T) 2 Aa(t)-r^Aa(t), (19) 

Acr(t) = a(t) - a cq (e), 

where cr oq (e) is the equilibrium value of the condensate which is non-zero for k = 0. The value of a cq is 
determined by minimizing the free energy, while the mass is related to the curvature of the free energy near the 
minimum, 



|-F(<t,T) = 0, ml(T) = -lf, 

ocr do 



—F(a,T) = 0, ml(T) = ^F(o,T). (20) 



One can see the equivalence of Eq. (|19fl with the differential equation for the harmonic oscillator after performing 
the following substitutions, 

fch.o./wh.o. — ► m 2 ., 7h.o./ m ii.o. — ► T, (21) 

where 7h.o. is the drag coefficient for the harmonic oscillator, fch.o. is the spring constant and mh.o. is the particle 
mass. The equivalence with the harmonic oscillator is shown in Appendix ^ along with a derivation of the 
same result from the perspective of linear response theory assuming a Langevin force added to the equations of 
motion. 

For the harmonic oscillator, the mean value of the position x is altered if the equilibrium position is moving. 
The amount of the change was consistent with the drag force jdxeq/dt being equal and opposite to the restoring 
force kbx. The corresponding result can be derived for the damped Klein-Gordon equation, 

7h.o. dx cq , T da eq 

Sx = ~Z aT> Sa= 2~FF\—]r> ( 22 ) 

fch.o. etc mtrx-L ) dt 

where 5a is the mean offset from the equilibrium value. Thus, m 2 a determines the restoring force, while T 
describes the drag. Finite-size effects could be estimated by replacing m 2 with m 2 + k 2 , where k 2 would be 
given by the finite size, k ~ 1/L. 

The resulting bulk viscosity is, 
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r d<J ca .„„. 

— -^s. (23) 



2 ds 



The bulk viscosity is then large for energy densities where a is rapidly varying, or for when m a is small, i.e., 
the critical region. 

III. BULK VISCOSITY IN THE LINEAR SIGMA MODEL 

Both of the last two sources of viscosity described in the previous section can be of special importance during the 
chiral transition. First, since masses are changing suddenly near T c , chemical abundances should easily stray from 
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FIG. 1: For the linear sigma model, the sigma field and mass are shown as a function of the temperature in the left panels. 
Near T c , the masses fall to zero and the mean value of the field changes rapidly, which gives rise to a sharp peak in the bulk 
viscosity. The pressure and temperature are displayed in the right-side panels as a function of energy density. The transition 
is sharper for g = 3.4 (solid lines) which is near to the critical value, than for g — 3.0 (dashed lines) which results in a smooth 
cross- over. 



equilibrium. Secondly, the mean field, which is zero above T c suddenly changes, and given the small masses in this 
region, large bulk viscosities are expected. 

As an example, we consider a simple example of a linear sigma model, where the coupling of the sigma field to the 
quarks provides the quark mass 



simple . 

Urn 



H 



1 



.VV+i^-^ + m^f 



h q a 



•ffquarkstVn = ffCr), 



(24) 



assuming only up and down flavored quarks. The resulting equation of state and values for m CT and a are displayed in 
Fig. \I\ior A 2 = 40. For couplings g < g c — 3.554, the transition is a smooth cross-over, while for g — g c the transition 
is second order, and for g > g c a first-order phase transition ensues with T c = 172 MeV. The values for g c and T c 
differ from 11] in that we neglect Fermi statistics here and also use a different value for A 2 . From Fig. ^ one can see 
that m a becomes small in the same region that the field rapidly changes, which results in a peak in the bulk viscosity. 
For g = g c , ra 2 a ~ \T — T c \, and the viscosity behaves ~ 1/|T — T c \. Since this is a mean field calculation we expect 
the critical exponent characterizing the singularity to differ from a more realistic quantum calculation. 
To calculate the bulk viscosity, we must evaluate the expressions from the previous section for 5n and 5a, 



5n 
5a 



dn oq 

-T c hem—j—SV ■ U, 

as 
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r da c , 
m 2 ds 
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These translated into changes of the pressure, which can be calculated as a function of a, T and the chemical potential 



rn dP r dP ^ <9P r 

SP =^ t+ 0f ST + l)a- Sa > 



(26) 



where the last term is zero since a is chosen to minimize the free energy. The change in the pressure can then be 
found by solving for S^jl and 5T using knowledge of Sn from Eq. I|25(l and the condition that e is fixed, 



r de de crTl de . 

6 *=8-^ + df ST+ 8a- 5 " = °> 

dn dn dn 

Sn= d; Sfl+ df ST+ d^ Sa - 



(27) 
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FIG. 2: The bulk viscosity, scaled by the energy density, is displayed in the upper panels for the linear sigma model. The peak 
at T c is due to the a field's inability to reach equilibrium, while the peak at low energy density is due to falling away from 
equilibrium. Viscous terms are larger and sharper for couplings close to the critical coupling (g c — 3.554). Lower Panel: For a 
Bjorken expansion (V ■ u = 1/r) the pressure is plotted alongside the Navier-Stokes expression, Tu — P — BX7 ■ u. Since the 
Navier-Stokes expression is only valid for small changes to the stress-energy tensor, the dynamics of the mean field should be 
handled explicitly if the corrections are large as are illustrated in the lower left-hand panel. 



Solving for dfi and 8T, 



<5/j = 
ST = 



(de/dT)5n+[(de/da)(dn/dT) - (de/dT)(dn/da)}Sa 
(de/dT)(dn/dn) - (de/d^)(dn/dT) ' 

(de/dn)8n+ [(de/dcr)(dn/dn) ~ (de/d^)(dn/da)]Sa 



(28) 



(de/dn)(dn/dT) - (de/dT){dn/dii) 
The partition function and pressure for the quarks can be calculated analytically for fixed a, 

TlnZ 24 



-Pquarks(m = g<T,T) 
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— {m 2 T 2 K (m/T) + 2mT 3 K 1 (m/T))} , 



2tt 2 



(29) 



where the factor of 24 is the number of degrees of freedom. Analytic expressions for the derivatives in Eq. I|28|) can 
also be found since derivatives of Bessel functions yield Bessel functions. 

Since equilibrium values are functions of T, the derivatives with respect to the entropy density in Eq.s l)25H can be 
expressed in terms of T, 



da, 



eq 



ds 

dn cq 

ds 



da cq /dT 
ds/dT ' 

dn e q/dT 
ds/dT 



(30) 



The derivatives of n cq and <7 e q in Eq.s l|30(l must be found numerically using the constraint that dP/da = at fixed T. 
After substituting the expressions for da cq /dT and dn cq /dT into Eq.s (|28|l . then substituting the resulting expressions 
for 5/j, and 5T into Eq. I|26|l . the resulting expression for 5P is linear in V • u, allowing the determination of the bulk 
viscosity using SP — — B\7 ■ u. 

The bulk viscosity was calculated assuming that the width T — 400 MeV, and that the chemical equilibration time 
scales inversely with the density and r c hem=l fm/c for a density of one quark per fm 3 . 

For a Bjorken expansion V • u = 1/r, and assuming an isentropic expansion starting with e = 8 GeV/fm 3 at r = 1 
fm/c, we calculated both P and B as a function of r. To illustrate the size of the effect, we display both P and 
the Navier-Stokes expression Tu = P — BV • u as a function of the energy density for a Bjorken expansion in Fig. 
121 The effect is certainly dramatic for g ~ g c , but since the Kubo relations are based on linear response theory, 



the larg e effect invalidates the underlying assumptions of Navier-Stokes hydrodynamics. We expect Israel-Stewart 
[l3L [Tj, HJ| [Til [13 equations for hydrodynamics to result in moderated effects compared to Navier-Stokes, though 
they should give identical results if the corrections are modest. If the effects are also large in Israel-Stewart solutions, 
one should consider treating the dynamics of the mean field explicitly along the lines of [llj , where the equations of 
motion for hydrodynamics and for the the mean field were solved in parallel. 

IV. SUMMARY 

The simplicity of the Kubo relations, Eq.s J7J) and Eq. JSJ, masks the wide variety of physical sources of viscosity. 
The one common aspect of the various sources is that non-zero equilibration times or non-zero interaction ranges can 
always be identified. We have focused on bulk viscosities associated with the chiral transition. In general, one would 
expect such effects whenever a system needs to rapidly rearrange its basic structure. In this sense these effects have 
much in common with super-cooling or hysteresis. In the case of a first order transition where super-cooling occurs, 
the departure from equilibrium is not proportional to the rate at which the system is changing, and the language of 
Navier-Stokes hydrodynamics is inappropriate. 

The peaking of the bulk viscosity near T c is in stark contrast to the behavior of the shear viscosity for many fluids 
which comes to a minimum near T c J18J . In |lq convincing physical arguments are presented that the shear viscosity 
for the deconfinement transition also comes to a minimum near T c . If the source of the viscosity is mainly due to the 
system's failure to equilibrate a scalar quantity such as the a field, one physically expects the singularity to be confined 
to the bulk viscosity. However, it is of interest that models of binary fluids also make predictions of a singularity in 
the shear viscosity near T c as described in |l9| . where physical arguments are made by thermodynamically linking 
the diverging correlation length to a divergence in viscous forces. One lesson from the study of critical phenomena is 
that critical exponents inferred from mean-field models like those discussed here will likely be incorrect, even though 
the qualitative behavior might be well reproduced. 

The implications for dynamics should be that the matter accelerates more quickly due to the higher gradients in 
T xx when the interior energy density is above the critical region. Once the matter flows into the viscous region of 
energy densities, there should be a slowing down and a reduction of surface emission. This trend would be in the 
right direction to explain HBT measurements which show a rapid expansion with a sudden disintegration |20j, but 
the potential magnitude of the effects are not yet known. 

Finally, we re-emphasize that if one were to solve for the evolution of the mean fields or chemistry alongside solving 
the hydrodynamic evolution equations, one could forego incorporating these effects through the viscosity coefficients. If 
the stress-energy tensor is strongly affected, the proper conclusion may be that rather than absorbing these effects into 
viscous hydrodynamics, one should treat non-equilibrated degrees of freedom explicitly. For instance, the dynamics of 
the a field can be calculated in parallel to the hydrodynamic equations of motion as was done in |ll| , though it would 
be important to incorporate damping effects which were neglected there. Furthermore, chemical non-equilibration 
might be accounted for by solving for the time evolution of chemical abundances, then expressing the pressure in 
terms of the resulting non-equilibrated densities. 

APPENDIX A: EQUIVALENCE BETWEEN DAMPED HARMONIC OSCILLATOR AND MEAN FIELD 

EQUATIONS 

The differential equations for the harmonic oscillator and for the Klein- Gordon equations, 

m-^ = - 7 ^-k(x-x (t)) + F(t), (Al) 



a 



da 



^ = -r^-™ 2 >-- eq ) + m 

are equivalent after the substitutions x ^ a, m <-> 1, 7 <-> T and k «-» m^. Here, F(t) is an external random driving 
force, i.e., a Langevin force. After solving the more physically intuitive harmonic oscillator, a simple substitution will 
provide the answer for the Klein-Gordon equation. 

Here, we consider the harmonic oscillator where the minimum of the potential moves with a velocity vq, xq = vot. 
Defining x' = x — x (t), one can find a differential equation, 

d x* dx f 

m—= 1 —-kx' + F(t)- 1 v Q . (A2) 



8 

Thus, the effect of moving the potential is a drag force 7^0, which causes the mean position to be offset from the 
center of the potential by an amount, 

5x' = -jvo/k. (A3) 

This illustrates the importance of the drag term, and the irrelevance of the mass term or the strength of the random 
term. 

The same result, Eq. (IA3|1 can be derived from the classical expression for linear response theory, 

(8x) = 13 fdt'(x(t = 0)dt>V(t')}, (A4) 



which leads to 



V(t) — ~k[x — xo(t)] 2 — -kx 2 rs —kxvot, 



,0 

(5x) = (kvo/T) / t'dt'(x(0)d t >x(t')) (A5) 

= -{kv Q /T) / dt'(x{0)x(t')). 
J —00 

The solutions can be written in terms of the external force in frequency space using Green's functions, 

x{lj) = (J){uj)F{uj), (A6) 

*M = 2 -■ -Z - (A7) 

The correlations in time can be written as: 

(af(0)a;(t)) = f dt' dt" (j>{-t')(f>{t - t")(F(t')F(t")), (A8) 

where for a Langevin force, (F(i)F(t')) — LS(t — t'). 

(x(0)x(t)) = L I dt'cf>(-t')<t>{t-t') (A9) 

= — [ dcje- iut (j)(u;)(j)(-uj). (A10) 

2n J 

The integration over time then is then 

r ° L 

dt(x(0)x(t))=L\^ = 0)\ 2 = ^. (Afl) 

The strength of the Langevin force L can be determined by calculating (x(0) 2 }, then using the equipartition theorem. 

(x 2 ) = L J dt'[^-t')] 2 = ±- J du\<f>{oj)\ 2 (A12) 

l r, r 1 1 

dio 



2nm 2 J \ (u) — i"f/2m — Q) (u> — i"f/2m + f2) 
1 1 



(w + ij/2m -Q) (lo + ij/2m + 0) 
L 

2fc^' 



, il = y/ik/m) 2 + (7/2m) 2 



which when combined with the equipartion theorem, (l/2)k(x 2 ) = (1/2)T, gives 

L = 2jT. (Af 3) 



As expected, the strength of the random force depends on the temperature and the damping, and is independent of 
the mass or spring constant. When combined with Eq. (|A4() yields 

(5x) = -^ (A14) 

which agrees with the simple statement above that on average the position of the mass lags behind the minimum 
of the potential by an amount such that the drag force, 7^0, cancels the restoring force of the spring, kSx. It is 
remarkable that the result is independent of either the mass or the temperature. 

The zero mass limit can be linked to the diffusion equation. To make this connection, we consider a random force 
F which acts for a time St. In the zero-mass limit of Eq. (|A1|) particles move by an amount, 

(k F\ 
Sx = -[-x + —)St, (A15) 

\7 7/ 

which means the distribution changes due to the underlying translation, 

«/W --/Ml* -fc^ + i*^ (MO, 

k utf \ , **_ .0/(3) , 1 (F 2 + k 2 x 2 )6t 2 d 2 f(x) 

= -Stf(x)-\ 5t— = h- 5 ^l - ' 

7 7 ox 2 7^ ax z 

where terms linear in F have been discarded assuming that the sign of the force is random. Taking the limit 6t — ► 0, 
and assuming the Force is of the Langevin form, F 2 St = L, one finds, 

o_m = k m+ te^ +D »m (A17) 

at 7 7 ox ox z 

L _T 

V ~~ 7' 



D = 



Thus, if the dynamics of a field are described by either a first-order or second-order differential equation, the problem 
can be mapped to the harmonic oscillator. Furthermore, if it is first order as is the case for chemical equilibration, it 
can be mapped to the diffusion equation. 

After making the substitutions described in Eq. IjAljl . the result translates to: 

The similarity between the equations of motion for a harmonic oscillator and for fields allow for an easier understanding 
of the terms in the Klein-Gordon equation. The field's effective mass, m 2 , plays the role of the spring constant and 
provides the restoring force pushing the field towards equilibrium, whereas the width T provides an effective drag 
force which impedes the field from maintaining equilibrium. This might be opposite to one's intuitive expectation: 
that a large width T would indicate a quick decay of any excitation of the field. 
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